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Abstract: Clusters of galaxies, such as the Virgo cluster, host enormous quantities of dark matter, making them prime 
targets for efforts in indirect dark matter detection via potential radiative signatures from annihilation of dark matter 
particles and subsequent radiative losses of annihilation products. However, a careful study of ubiquitous astrophysical 
backgrounds is mandatory to single out potential evidence for dark matter annihilation. Here, we construct a multi- 
wavelength spectral energy distribution for the central radio galaxy in the Virgo cluster, M87, using a state-of-the-art 
numerical Synchrotron Self Compton approach. Fitting recent Chandra, Fermi-LAT and Cherenkov observations, we 
probe different dark matter annihilation scenarios including a full treatment of the inverse-Compton losses from electrons 
and positrons produced in the annihilation. It is shown that such a template can substantially improve upon existing dark 
matter detection limits. 

Keywords: dark matter, clusters of galaxies, high-energy emission, extragalactic jets, M87 



1 Introduction 

The presence of dark matter in the universe can be astro- 
physically observed through the spectrum of acoustic reso- 
nances in the primordial plasma imprinted on the spectrum 
of anisotropics in the cosmic microwave background radi- 
ation, and its gravitational effects on visible matter. This 
dark matter could be weakly interacting massive particles 
(WIMPs) with masses at the electroweak symmetry break- 
ing scale ^ O (^{GfV2}^^^ = 246 GeV, freezing- 
out in the early universe. Their reUc density would be in 
agreement with the observed ridmh'^ = 0.1120 fl]. This 
is termed the "WIMP Miracle". Moreover, stable neutral 
particles at the electroweak symmetry breaking scale are 
predicted by supersymmetric extensions of the Standard 
Model. 

In the present day universe, annihilation interactions in 
regions with high mass densities still occur and can 
lead to the production of energetic neutrinos, electrons, 
positrons and gamma rays which inverse-Compton scatter 
with starlight photon fields or the cosmic microwave back- 
ground resulting in signals potentially detectable across the 
soft to hard X-ray energy band. Limits on this dark matter 
annihilation (DMA) emission have been studied thus far 
in the Galactic Center |2|, extragalactic gamma ray back- 
ground lis, clusters of galaxies ID and by focusing on 
dwarf galaxies in the Milky Way halo |I5]|6|. Searching 
for signatures of secondary gamma rays and decay prod- 



ucts lends insight into the thermally averaged annihilation 
cross section and mass of the dark matter particle. 
The first radio galaxy detected in the TeV energy range, 
M87 is one of the best studied in its source class with a 
jet resolved by optical, radio and X-ray observations. This 
cD galaxy is the nearest radio galaxy emitting very high 
energy (VHE) gamma rays and is known for its proxim- 
ity to the Earth (D — 16 Mpc) and exceptionally bright 
arcsecond scale jet Q. M87 is an excellent laboratory in 
which to study dark matter because of its location in the 
central high mass density region of the Virgo cluster. In this 
work, the imprint of WIMP annihilation radiation from the 
halo structure of M87 in the form of prompt and inverse- 
Compton gamma rays is investigated, and dark matter par- 
ticles are constrained, via a multiwavelength analysis of the 
spectral energy distribution of M87. 

2 Observations and Analysis 

Observations from the Chandra X-ray Observatory, Fermi 
Large Area Telescope (Fermi-LAT), Major Atmospheric 
Gamma-ray Imaging Cherenkov Telescope (MAGIC), and 
High Energy Stereoscopic System (H.E.S.S.) Cherenkov 
array are used in this analysis. In contrast to the model 
presented by Abdo et al. Q, historical data in the low- 
energy range from the radio-to-optical regime which can be 
contaminated by dust and starlight and are thought to orig- 
inate from different regimes in the source than the VHE- 
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Figure 1: (a) Spectral energy distributions of M87 using the Synchrotron Self Compton model (light grey line); and with 
inclusion of dark matter annihilation model (dark line). The SSC model alone produces a value of 27.5 whereas an 
improvement to 14.4 is achieved using the dark matter annihilation model for a particle of mass 4.7 TeV and a boost 
factor of 812. The data points used in this analysis are shown, with inverse-Compton and prompt pion gamma peaks also 
illustrated. There is a noticeable gap between observed data and the SSC model in the MAGIC regime which is closed 
by the SSC+DMA model, (b) Island plot showing result of analysis for the combined M87 fit. The boost factor / is 
introduced as a free parameter to account for the enhancement due to sub-halo clumping; m denotes the generic WIMP 
mass. 



emission are not included in the fit but rather used as bench- 
marks that the emission should not exceed. Although non- 
contemporaneous, the data set represents a long-term av- 
erage by excluding significant flaring events and thus rep- 
resents the steady low state spectral energy distribution of 
high-energy emission in M87. 

The Synchrotron Self Compton (SSC) process occurs when 
electrons accelerated in the magnetized plasma jet emerg- 
ing from the ergosphere of a central black hole reach ultra- 
relativistic speeds and emit synchrotron photons which 
then inverse-Compton scatter off the original source pop- 
ulation of electrons. This generates a high-energy spec- 
trum that extends into the TeV regime. The Synchrotron 
Self Compton model of Riiger et al. [,8J is applied here 
to the high-energy data with a focus near the upper cut-off 
of the spectrum so that the Compton scattering cross sec- 
tion is treated appropriately in the Klein-Nishina regime. 
This model has relativistic electrons injected into a spher- 
ical emission region with a randomly oriented magnetic 
field, moving up the jet with relativistic speed. Taking a 
bulk Lorentz factor F = 2.3 (Doppler factor 5 = 3.9) as 
in [71 we vary other parameters so the maximum Lorentz 
factor Fniax = 10*, magnetic field S = 3 G, normalization 



factor K = 10 cm ''s 



and the radius of the emitting 



source R = 3.5 x 10^^ cm to obtain a differential slope 
of s = 2.2 for the injected electron power law distribution 
with exponential cut-off. 

A dark matter model is then added to this SSC model. 
When dark matter neutralinos annihilate, they produce 
heavy quarks, leptons and W bosons. During their subse- 
quent hadronization, they decay mainly into pions. Prompt 
pion emission is the emittance of very high energy gamma 



rays from the decay of neutral pions. The charged pions de- 
cay into electrons and positrons which can up-scatter cos- 
mic microwave background photons through the inverse- 
Compton mechanism also to very high energies. This 
model is referred to here as SSCh-DMA. 
The differential flux of gamma rays produced by the de- 
cay of annihilation products in M87 is given by Equation 
(1). {(Jav) is the thermally averaged annihilation cross sec- 
tion, the mass of the dark matter particle, and / the 
boost factor that accounts for enhancement due to sub-halo 
clumping of dark matter in the M87 halo. D is the distance 
from the Earth to M87, p{r) the dark matter density profile 
and the gamma photon spectrum coming from the de- 
cay of neutral pions from hadronization in the annihilation 
process (prompt pion emission). 
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The differential flux of high-energy photons produced by 
inverse-Compton scattering of charged particles off the 
cosmic microwave background is given by Equation (2). 
P {E, E') is the differential power emitted into photons 
of energy E by an electron or positron with energy E' . 
b [E') is the total rate of electron/positron energy loss due 
to inverse-Compton scattering as in |9|, and is the 
sjjectrum of secondary electrons and positrons with energy 
E. 
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Finally a Navarro-Frenk- White density profile is assumed 
for the dark matter halo so that p{r) = pNFw{r) ifTOl . a 
generic species of annihilating WIMPs with rest mass in 
the GeV - TeV range, and a thermally averaged annihila- 
tion cross section {aAv) = 3 x lO^'^'* cm'^s^^. 10^ real- 
izations of generic WIMP neutralinos motivated by Super- 
symmetry were generated to obtain the average secondary 
electron/positron spectrum The gamma ray emission 
spectrum due to the decay of neutral pions in the annihi- 
lation process was generated numerically using the Dark- 
SUSY code ifTTl . The prompt pion emission and inverse- 
Compton mechanism reflect the dark matter halo distribu- 
tion while usual Synchrotron Self Compton processes orig- 
inate from the jet. 

The best-fit combined model of AGN emission and the 
complete dark matter related emission is determined using 
a x^-test. See Figure 1 . 



3 Results 

In addition to closing the noticeable gap between the ob- 
servational data and the SSC model in the high-energy 
regime, introducing a WIMP annihilation component to 
the Synchrotron Self Compton emission model of multi- 
wavelength observations of M87 produces a minimal 
value and thus an improved agreement with the data set. 
The best-fit employed a neutralino mass — 4.7 TcV 
and boost factor / = 812 both within the reach of numeri- 
cal models and modern collider experiments. 
Assuming annihilation interactions of WIMPS with TeV- 
scale masses in the dark matter halo of M87 are responsi- 
ble for prompt gamma ray and inverse-Compton emission 
components, this method can be used to constrain viable 
dark matter models as new observations are being made 
eg. through Fermi-LAT and MAGIC telescope observa- 
tions. Furthermore, insight into the discrimination between 
diffuse halo emission and jet emission as well as sub-halo 
clumping effects can be found through ongoing observa- 
tional and theoretical investigation of the spatial distribu- 
tion of the gamma ray emission in M87. The gamma ray lu- 
minosity due to dark matter annihilation in M87 can reach 
the same level as the emission from its relativistic jet for 
reasonable boost factors. 

Since the high-energy emission of M87 shows blazar-like 
variability |12|, the constraints for dark matter models 
can be improved significantly with future contemporane- 
ous multiwavelength data, measuring the spectral energy 
distribution during low flux states. Whereas the gamma 



rays from the jet show correlated broad-band variability, 
the quiescent components produced by pions and inverse- 
Compton scattering electrons should remain steady. The 
DMA prompt pion and inverse-Compton emission result in 
a characteristic double-hump structure superimposed on 
the SSC emission which could in future be used as a distin- 
guishing characteristic. A dark matter induced component 
to the spectral energy distribution would stand out as steady 
emission. Thus there is potential to recognize this during 
low state emission of the active galactic nucleus and possi- 
bly extract dark matter properties such as the WIMP mass, 
annihilation cross section and boost factor due to substruc- 
ture. Among the steady background components, the only 
other potentially confusing source of gamma rays is due 
to cosmic rays traveling through the intracluster medium. 
This component should, however, be well distinguishable 
owing to its smooth power law and different spatial profile 
lfT3l . Although the significance of this result is not suffi- 
cient to claim evidence of a specific particle it is encour- 
aging further studies, especially since previous deep expo- 
sures of other targets like the Galactic Center 1 14| have also 
provided tentative hints of Te V-scale WIMP annihilation. 
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